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ABSTRACT 

We have observed the supernova remnant MSH 15—52 (G320.4— 1.2), which contains the gamma-ray 
pulsar PSR B1509— 58, using the CANGAROO-III imaging atmospheric Cherenkov telescope array 
from April to June in 2006. We detected gamma rays above 810 GeV at the 7 sigma level during 
a total effective exposure of 48.4 hours. We obtained a differential gamma-ray flux at 2.35 TeV of 
(7.9± 1.5stat ± l-7sys) X 10-13 cm-2 s'^ TeV'^ with a photon index of 2.21 ±0.39stat ±0.40sys, which is 
compatible with that of the H.E.S.S. observation in 2004. The morphology shows extended emission 
compared to our Point Spread Function. We consider the plausible origin of the high energy emission 
based on a multi-wavelength spectral analysis and energetics arguments. 

Subject headings: gamma rays: observations - ISM: individual(MSH 15—52, G320.4— 1.2) - supernova 
remnants - pulsars: individual(PSR B1509-58) 



1. INTRODUCTION 

Recent imaging atmospheric Cherenkov telescopes 
(lACTs) have achieved remarkably high sensitivity in 
the very high energy gamma-ray band. This is well 
illustrated by the Galactic pla ne survey carried out 
by the H.R S.S. coUaboration (jAharonian et al.l l2006l : 
iHoppd 120071 ). Many of the discovered TeV sources 
are associated with puls ar wind nebulae (PWNe) (e.g., 
iGaensler fc Sla^[20M ). which are now established as 
the most populous category among G alactic TeV sou rces: 
18 PWNe have been found so far (|Hintonl |2007[ ) and 
a portion of 21 unidentified galactic TeV sources could 
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be PW Ne as well (jAharonian et all l2007bl : iFunk et all 
l2007bD . 

The Crab nebula, the prototype PWN, has been 
observed in every accessible wave band. The nebula con- 
tains the Crab pulsar, which has the highe st spindown 
energ y loss of known gamma-ray pulsars (jThompson I 
l2003l ). However, the radiation mechanism remains 
controversial for the Crab nebula and for other PWNe: 
it has been argued that the Crab nebula's broadband 
spectral energy distribution (SED) can be well explained 
by synchrotron emission in an average magnetic field of 

0.1 — 0.3 mC for radio to soft gamma-ray bands, and 
by inverse Compton (IC) scattering of synchrotron, IR, 
millimeter and cosmic microwave background (CMB) 
photons for the hard ga mma- ray band up to 100 TeV 
("Aha ronian k, A tovan' '1 996t lAtovan fc AharonianI 
1996: lAharonian. Atovan fc Kifund 119971 : 

Aharonian fc AtovanI 19981 : lAharonian et al.l l2004j) . 



On the other hand, it was claimed that gamma rays 
from the decay of neutral pious, produced by hadrons, 
contribute significantly at ^ 10 TeV for the Crab 
nebula (Bednarek & B artosikll2003f ). Consideration has 
also been given to whether the energy source of TeV 
gamma-ray emission from an SNR containing a pulsar 
(composite SNR) is the pulsar's spindown energy or 
related to the supernova explosion (e.g. iFunk etldl 
l2007b), and the efficiency of energy conversion to the 
particle acceleration for such models. 

PSR B1509— 58 has the third highest spindow n en- 
ergy loss after t he Crab pulsar (iThompsonl l2003f ) and 
PSR J1833-1034 (jCamilo et al.ll2006f l in th"eGalaxv. and 
its nebula has also been well studied across the electro- 
magnetic spectrum. TeV gamma-ray observations, com- 
bined with those at other energy bands, provide addi- 
tional information which may lead to solutions for the 
above problems and a unified comprehension of pulsar 
and nebula systems. PSR B1509— 58 was detected in 
the radio supernova remnant MSH 15-52 (G320.4-1.2) 
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(| Caswell et alj [l98ll ). initially as a 150-ms X-ray pu l- 
sar with the Einstein satellite (jSeward fc HardenlflQSl ). 
which was confi rmed by later X- ray /soft gamma- 
ray observations (iKawai et al. 1993c iMatz et al. 1994; 
Q\xmi ct al. 1994; Saito ct al. 1997; Mars den et al.lll99'^ : 
Cusumano et al. 2001; For ot ct al. 2006,) . Subsequently 
its pulse period was detected at radio frequencies 
(jManchester et al]|1982f). and at soft gam ma-ray ener- 
gies with COMPTEL (|Kuiper et all [19991 ). The pulsar 
was detected above 30 MeV by EGRET at the 4.4cr level 
with some suggestive, if not statistically c ompelling, ev- 
idenc e of modulation at t he pulsar period (iKuiper et al.l 
I1999D. Although optical (IShearer et al.l ll998^ and near- 
IR (jKaplan fc Moonll2006D searches found possible pul- 
sar counterparts, the pulse period was not detected. As 
one of the most energetic young pulsars, PSR B1509— 58, 
has been particularly well studied at radio wavelengths. 
A detailed timing analysis yielded a period derivative 
of P = 1.5 X 10~^^, and a high spin-down luminos- 
ity of £^ = 1.8 X 10'^^/45 ergss"'^, where £45 is the mo- 
ment of inertia in units of lO^^g cm~^. A braking in- 
dex of n = 2.84 was measured, corresponding to an 
age, assuming an initial period of Pq <C P, of r = 
(P/(n - 1)P)[1 - (Po/P)"-^] ~ (P/(n - 1)P) - 1700yr. 
A large dipole surf ace magnetic field of i? = 1.5 x 10^^ G 
was also inferred (jKaspi et al.l 119941 : iLivingstone et al.l 
I2OO5I ). However, in contrast to other young pulsars, no 
pulsar glitch has been observed to date. Ini tially, the age 
of the SNR was esti mated to be 6 — 20 kyr (ISeward et al.l 
Il983| ) or - lOkyr (jVan den Bergh fc Kramper I Il984( l 
promptin g debate about the dis agreement with the pul- 
sar's a ge. iBlandford fc Romanil |l988) and lGvaramadze I 
(|2001h . for instance, suggested the pulsar age was actu- 
ally > 20 kyr. The latter was based on the assumption 
that pulsar's braking torque was enhanced by the inter- 
action between the pulsar's magnetosphere and circum- 
stellar dense clumps. 

The radio morphology of MSH 15—52 consists of south- 
east and northwest shells. The latter, ~ 10' from the 
pulsar, spatiall y coincides with the Ha nebula RCW89 
(|Rodgers et"all Il960f ). and iGaensler et"aI1 (|1999( l con- 
cluded that MSH 15-52, PSR B1509-58 and RCW89 
were associated systems. The distance, derived from an 
HI absorption measurement, is 5.2 ± 1.4 kpc, consistent 
with the value of 5.9 ib0.6kpc determined fr om the pulsar 
dispersion measure ()Tavlor fc Cordes |[l993i) . We adopt 
d = 5.2 kpc throughout this paper. 

Symmetric Jets, sini ilar to the Crab pulsar 
(Brinkmann ct a l. | | 1985D and the Vela pulsar 
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were observed by ASCA 
EOS A T (iTrussoni et al.l [l99fih 



and Chandra (|Gaensler et al.l l2002t l. The jet directed 
towards the northwest was observed to terminate at 
RCW89. Precise Chandra observations revealed the 
seq uential heating of t he knots in RCW89 by the pulsar 
jet (|Yatsu et al.l l2005i ) . The high resolution Chandra 
image also showed the arc structure where the pulsar 
wind ma y be terminate d , and the diffuse pulsar wind 
nebula (jCaensler et al.l l2002| ). which emitted non- 
thermal synchrotron radiation. However the PWN has 
not been observed at other wavelengths such as IR or 
optical. Although IRAS found an infrared source, IRAS 
15099-5856, spatially coincident with PSR B1509-58, 



the IR emissi on was therm al and therefore not related 
to the P WN (lArendtlll99lD. A faint radio structure was 
detected (jCaensler et al.lll999l) b ut its flux density wa s 
obtained only within a large error (jCaensler et al.ll2002f) . 
A larger extent than that of X-ray was expected due 
to the difference of the cooling lifetime so that it was 
thought to be partially hidden by the bright RCW89. 

The Ginga LAC discovered single power-law emis- 
sion up to 20keV with a photon index of 2, indi- 
cating synchrotro n emission and the existence of accel- 
erated electrons (jAsaoka fc Kovanial Il990t) . The syn- 
chrotron nebuk^_S2ectrijin^ also mea sured by EIN- 
STEI N (ISeward et al.l[l983h. EX OSAT (jTrussoni et all 
[1993) and RXTE (iMarsden et aLl[l99l . Bepyo SAX de- 
lected nonthermal emission from 1 keV u p to 200 keV 
with a photon index of F = 2.08 ± 0.01 (jMineo et all 
120011 ). while recent observations with INTEGRAL IBIS 
found a possible (2.9g) spectral cut off at ~ 160 keV 
(jForot et al.ll2006l ). Since high energy electrons were 
demonstrated to exist, very high energy gamma-ray emis- 
sion was predicted via IC scattering w ith CMB photons 
(jdu Plessis et al.l[l995l : [Harding I [l996h . CANCAROO- 
I also suggested a possible VHE gamma-ray detection 
of ^ 10% of the Crab flux above 1.9 TeV, a ssuming a 
spectral photon index of 2.5 (jSako et al.ll2000D . H.E.S.S. 
subsequently reported extended VHE gamma-ray emis- 
sion along with the pulsar jet. Their TeV morphol- 
ogy showed a good coincidence with X-ray images, in- 
dicating that the TeV gamma-rays originate from the 
inverse Compton scattering of relativistic electrons. It 
was pointed out that IC scattering of the CMB could 
not account for most of the TeV gamma-ray flux with 
an assumed magnetic field of 17 /iC, which indicated the 
contribu tion of IR photons as seed photons for t he IC 
process (jAharonian et al.|[200l [Khelifi et al.|[2005l ). The 
necessity for target photons in addition to the CMB has 
recently been sugg ested for other Galactic sources also 
(jKhelifi et al.ll2005l : lHinton fc Aharonianll2007f ) . Here we 
report the results of TeV gamma-ray observations with 
the CANGAROO-III telescopes and consider the origin 
of the TeV gamma-ray emission based on a discussion of 
the energetics. 

2. CANGAROO-III OBSERVATIONS 

CANGAROO-III is an array of four lACTs, located 
at Woomera, South Australia (136°47'E, 31°06'S, 160m 
a.s.l.). Each telescope has a 10m diameter reflector 
which consists of 114 segmented FRP spherical r nirror s 
mounted on a parabolic frame (jKawachi et al.l |2001| ). 
The telescopes are situated at the c orners of a diamond 
with ~100m sides (Eiiomo t'oet al.l [2002). The oldest 
telescope, Tl, which was the CANGAROO-II telescope, 
was not used due to its smaller FOV and higher energy 
threshold. The imaging camera systems on the other 
three telescopes (T2, T3 and T4) are identical, w ith 427 
PMTs and a FOV of 4.0° (jKabuki et all [2001 . The 
PMT signa ls were recorded b y charge ADCs and multi- 
hit TDCs (jKubo et al.ll2003f l. The observations were 
made from April to June in 2006. The tracking posi- 
tions were offset by ±0.5° from PSR B1509-58 in dec- 
lination or in right ascension, and changed every twenty 
minutes, in order to suppress position-dependent effects 
on the camera due to bright (4.1 and 4.5 magnitude) 
stars. To trigger data recording, an individual telescope 
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was required to have more than four pixels registering 
over 7.6 photoelectrons within lOOnsec (local trigger), 
with a global trigger system then determ ining the coinci- 
dence of any two of the three telescopes (iNis hiiima et al.l 
l2005f ). We rejected data taken in bad weather conditions 
in which the shower event rate was less than 5 Hz or at 
zenith angles larger than 35°. Finally, the selected data 
were taken at a mean zenith angle of 30.1° with a cor- 
responding Point Spread Function (PSF) of 0.23° (68% 
containment radius). A typical trigger rate of 3-fold co- 
incidence is 12 Hz. The effective exposure time amounts 
to 48.4 hours. 

3. DATA REDUCTION AND ANALYSIS 

The basic analy sis procedures are described in detail 
in lEnomoto et al.l (2006) and iKabuki et al.l (|2007( ). Us- 
ing calibration data taken daily with LEDs, the recorded 
charges of each pixel in the camera were converted to 
photoelectrons. At this step we found 8 bad pixels out 
of 427 X 3 due to their higher or lower ADC conversion 
factors in these observations. These bad pixels were re- 
moved from this analysis period, which was also reflected 
in the Monte Carlo simulations. Then every shower im- 
age was cleaned through the following CANGAROO-HI 
standard criteria. Only pixels which received >5.0 pho- 
toelectrons were used as "hit pixels" . Then five or more 
adjacent hit pixels, with arrival times within 30nsec from 
the average hit timing of all pixels, were recognized as a 
shower cluster. 

We carefully studied the effect of the bright stars by 
monitoring the mean ADC counts of the pixels for which 
the stars entered the FOV. When the stars were within 
the pixel's FOV, the PMT hit rate increased significantly. 
At such times the "hits" for that pixel were dominated 
by the starlight triggers, which were dimmer than "hits" 
by air shower Cherenkov photons. The average ADC 
counts of the affected pixels were clearly reduced during 
such times. After the image cleaning procedure above, 
the averaged ADC counts became stable within the usual 
values over the whole run, indicating the effects of the 
stars were completely removed. 

Before calculating image moments — the "Hillas pa- 
rameters" (Hillas 198 5j) -— we a pplied the "edge cut" 
described in iKabuki et al.l ()2007[ ). If the brightest 15 
pixels in the image (or all pixels for images with less 
than 15 hits) were not located in the outer edge layer of 
the camera, we retained the event. The orientation an- 
gles were determined by minimizing the sum of squared 
widths with a constraint given by the distance predicted 
by the Monte Carlo simulations. 

Theii we a p plied the Fishe r Disc riminant method 
(jFished 119361 : lEnomoto et al.l l2006f l with a multi- 
parameter set of P = {W2,W3,W4, L2, L3, L4), where 
W and L are the energy corrected width and length, and 
suffixes identify the telescope. The Fisher Discriminant 
(FD) is defined as FD = a ■ P where c5 is a set of coeffi- 
cients mathematically determined to maximize the sep- 
aration between FD for gamma-rays and hadrons. 

For the background study we selected a ring region 
around the target, 0.2deg^< 9^ < 0.5 deg^, and obtained 
FD distributions for background, Ff,, and Monte Carlo 
gamma-rays, Fg. Finally we could fit the FD distribu- 
tions of the events from the target with a linear combina- 
tion of these two components. The observed FD distri- 
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Fig. 1. — FD distribution. The open circles show the FD ob- 
tained from the ON source region, 9^ < 0.1 deg^. The broken and 
solid histogram are the background and gamma-ray component es- 
timated by the fit procedure described in the text. The filled circles 
are the subtraction of the background from the ON source region. 
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Fig. 2. — The 8^ plot. Here, deg corresponds to the best-fit 
position of the centroid of the emission obtained by this work (see 
text). The hatched histogram represents our PSF derived from the 
Monte-Carlo simulation. 

butions, F, should be represented as F = aFg + {l — a)Ft 
where a is the ratio of gamma-ray events to total events. 
Here only a is optimized and the obtained FD distri- 
butions are shown in Fig. [TJ This analysis method was 
checked by an analysis of the Crab nebula data taken in 
December 2005. 

The reflectivities of each telescope, which are used in 
the simulations, are monitored every month by a muon 
ring analysis of calibration run individually taken by each 
telescope. We obtained relative light collecting efficien- 
cies with respect to the original mirror production time 
of 0.61, 0.64, and 0.67 for T2, T3 and T4, respectively. 

4. RESULTS 

The obtained 0'^ plot is shown in Fig.|2|with the PSF of 
our telescopes. Above 810 GeV we detected 427±63 ex- 
cess events under the assumption the TeV source was a 
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Fig. 3. — Morphology of gamma-ray like events, smoothed with 
a gaussian of a = 0.23°. Our PSF is also shown by a dashed 
circle(68% containment radius). The squares and the cross repre- 
sent tracking positions and the pulsar position, respectively. Solid 
contours show VHE gamma-ray emissi on as seen by H.E.S. S. and 
dotted contours by ROSAT0.6-2.lkeV l|Trussoni et al.ll996I V The 
region between thin dotted circles are used for the background 
study. 

point source {9^ < 0.06 deg^) and 582±77 events (7.6cr) 
within 6^ < 0.1 deg^, which corresponds to the size of 
the SNR. The TeV gamma-ray emission is extended, and 
the morphology of gamma-ray-hke events, smoothed by 
a gaussian with a = 0.23 deg, is shown in Fig.[3l The 
number of events were individuaUy estimated by the FD- 
fitting method in each 0.2° x 0.2° sky bin. When we 
evaluate the outer regions {9^ > 0.6 deg^), we must con- 
sider the gradual deformations of the FD distributions 
at larger angular distances from the target. Therefore 
we selected a ring with radii 0.2° < r < 0.4° centered 
on the evaluated region as the background. For the in- 
ner regions, 9^ < 0.6 deg^, the same background as used 
in Fig. [2] was adopted. The intrinsic extent of the TeV 
gamma-ray emission was estimated by the 2D Gaussian 
fit on our unsmoothed excess map. The inclination of 
the major axis is 61.3 ± 1.9°, measured to the north 
from west, which is compatible with the jet direction 
determined by Chandra, H.E.S.S. and INTEGRAL. The 
intrinsic source sizes along the major and minor axes 
were calculated to be 0.07 ± 0.07° and 0.21 ± 0.08°, re- 
spectively. The center of gravity of the TeV gamma- 
ray emission is at (R.A., Dec.)=(228°.486, -59°.235), 
which corresponds to the offset from the pulsar in (R.A., 
Dec.) = (0°.0030 ±0°.0076, 0°.10±0°.012). No signif- 
icant offset from the pulsar is then observed given our 
PSF, while the H.E.S.S. measurement showed the off- 
set of (0°.048,0°.022) in (R.A., Dec.) at the 3a level. 
Fig. m represents a reconstructed VHE gamma-ray spec- 
trum compatible with a single power-law: the value at 
2.35 TcV is (7.9±1.5stat±1.7sys) x 10"" cm-^s-^TeV"! 
with a photon index of 2.21 ± 0.39stat ± 0.40sys. The 
flux was measured within 9^ < 0.1 deg^. The rele- 
vant systematic errors are due to the atmospheric trans- 
parency, NSB fluctuations, uniformity of camera pixels, 
and light collecting efficiencies. In addition, the signal 
integrating region was changed from 9^ < 0.1 deg^ to 
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Fig. 4. — Differential flux of the whole nebula. Squares and 
triangles show the CANGAROO-III and the H.E.S.S. data points, 
respectively. The best fit power-laws are also shown by the solid 
and dashed line from this work and from H.E.S.S., respectively. 

0.2 deg^ in Fig.O which was included in the systematic 
error. The TeV gamma-ray extent and flux obtained 
by CANGAROO-III are consistent with those obtained 
with H.E.S.S. Our result indicates the TeV gamma-ray 
emission does not vary significantly between the H.E.S.S. 
observations in 2004 and ours in 2006, which is also 
consistent with the steady X-ray emission from the dif - 
fuse PWN over several decades (|Delanev et al.l l2006| ) . 
ISako et all (|2000f ) reported a marginal detection (4.1cr) 
above 1.9TeV of (2.9±0.7stat) x IQ-^^ ergscm^ s^i. For 
the sake of comparison, we integrated the differential flux 
of this work over the sa me energy band wi th a photon in- 
dex of 2.5 as inferred bv lSako et al.l (|2000f ). and obtained 
(1.6 ± 0.4stat) X 10~^^ ergscm"^ s~^, marginally outside 
the Icr error range of the CANGAROO-I result. 



5. DISCUSSION 

We discuss the possible origin of the TeV gamma-ray 
emission from the PWN for both proposed cases of its 
age: r = 1700yr an d 20kyr (jBlandford fc Romanil[l988l : 
lGvaramadze1l2001[ ). 

5.1. Upper limit on the global energetics 

In general, two alternatives of the energy source of a 
composite SNR could be considered, namely the pulsar's 
spindown energy and the supernova explosion (or the 
SNR). In the latter case, the energy could be as high as 
~ 10 ergs, or even up to a few times lOffergs for some 
SNelb/c and SNelln (|Nomoto et all 1200 1[ ). However, 
in the case of MSH 15—52, the extent of the intrinsic 
TeV gamma-ray morphology shows a good coincidence 
not with the shell of MSH 15—52 but with the jet and 
the PWN. Therefore we rejected the scenario in which 
the SNR is at the origin of the high energy particles. 
Hereafter we assume the pulsar's spindown energy is the 
global energy source. 

The time evolution of the pulsar period P{t) is de- 
scribed as 

P(i)=Po(l + -)'^, (1) 

V To/ 

where tq is a parameter called the initial spindown time 
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TABLE 1 

TO DEPENDENCE OF PARAMETERS. 



PSR 



TO [yr] Etot [ergs] Pq [msec] 



B1509-58 


30 


7.5 


X 


10« 


16 


(r=1700 yr) 


100 


2.1 


X 


10« 


31 




300 


7.2 


X 


1048 


53 




500 


4.6 


X 


1048 


67 




700 


3.6 


X 


10^8 


77 




1000 


2.8 


X 


1048 


87 


Crab 


30=^ 


1.1 


X 


1051 


3.3 




700'' 


3.7 


X 


1049 


19 



^ lAtovanl 1)19991 ) argued that the observed radio 
spectrum suggested ro = 30 yr. 
^ derived from eq. ijjjl. 

assuming both n and k in the braking equation 
-fcil" are constant (e.g.. iGaensler fc Slane I l2006f ) . 

The pulsar's spindown energy E{t) is then calculated as 

E{t)=EJl + -Y", (2) 



scale 



1 



where Eq is the initial spindown luminosity. By integrat- 
ing this, we obtain the total energy which the pulsar has 
lost over its age, r. 



Etnt — 



E{t)dt = 



EqTo 



l-a 

1 + -) -1 
To 



(3) 



This formula provides an upper limit to the energy 
supplied by the pulsar to the PWN as a function of 
the unknown parameter tq. For the sake of simplicity, 
energy loss through adiabatic expansion was neglected 
here. The dependence of E^tot on tq is listed in Table [TJ 
where the pulsar's moment of inertia was assumed to be 
10'*^gcm~^. The value of tq also determines the initial 
spin period Pq, using the current period P{t), as listed in 
Table [TJ For comparison, the Etot and Pq for the Crab 
pulsar were also listed in Table [H with P = 33 msec, 
P = 4.2 X 10-13 (iTavlor fc Co"rdes1 fl99l) , r = 950yr , 
E{t) = 5 X lO^^ergss-i and n = 2.5 (|Lvne et al.lll988D . 
The ii^tot of the Crab is about 10 times larger than that 
of PSR B1509-58 for the same tq. 
For T = 1700 yr, we can calculate tq as 



To 



P{r) 



(n-l)P(r) 



30yrs. 



(4) 



However, in the r = 20kyr case, these formulae are not 
applicable due to th e possibility of a time-dependence for 
k (Bland ford fc Rom ani 1988). It is hard to estimate an 
accurate i^tot; however, a larger Ptot than that estimated 
for T = 1700 yrs would be expected. 



5.2. Hadronic 



scenario 



Hadronic gamma-ray productio n in PWNe has been 
sugge s ted for the Crab n ebula (|Bednarek fc BartosikI 
12003"; 'Amato et al.' '2003") an d the Vela X region 
iBcdnarck 2007; Horns 2006; H orns et al.ll2007[ ). First 
we examine a neutral-pion decay model for the origin of 
the TeV gamma-ray emission. Fig. [5] shows the spectral 
energy distribution (SED). We assumed the population 
of accelerated protons to be expressed by a single power- 
law with an exponential cutoff, that is, dNp/dEp oc 

Ep^"" exp{—Ep/E^s,y^). We used only our data and the 
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A H.E.S.S. 
A CANGAROO-I 
* EGRET 
(> COMPTEL 
— 7t° decay 

pulsed component 




10"^ 10'- 10" 
Energy [eV] 

Fig. 5. — Spectral energy distribution and the model curve of 
neutral pion decay. The solid line represents the best fit model 
curve. COMPTEL and EGRET flux points were omitted from the 
fit since they contained pulsed emission from the pulsar. A dashed 
line cited from Cusumano et al. (2001) represents a model fit curve 
of the pulsed emission obtained hyBeppoSAX and COMPTEL. 



H.E.S.S. data llAharonian et al.ll2OO50 for the model fit 
since the COIVIPTEL and EGR ET data were the sum of 
pulsed and unpulsed emission (|Kuiper et al.lll999[) which 
may contain emission from the pulsar. The contribution 
of the pulsed emission is indicated by the dashed line in 
Fig. [51 derived from a fit to the BeppoSAX and COIVIP- 
TEL data of the pulsed emission (jCusumano 
(see also Fig. |S]) . Since these fluxes were apparently dom- 
inated by the pulsed emissio i i, we d idn't use them in the 
model fit. lAharonian et al.l ()2007aD put upper limits to 
the TeV gamma-ray pulsed emission of approximately 
one order magnitude below the total fiux, and this com- 
ponent was excluded from the fit. The best-fit curve is 
shown in Fig. [5] by the solid line. The power-law index 
and the cutoff energy were obtained to be 7p = 2.16±0.05 
and i?niax = 530 ± 399 TeV, respectively. The total en- 
ergy of high energy protons above 1 GeV, Wp, was cal- 
culated to be Wp = 3.2 X 10^i(n/lcm-3)-i ergs, which 
means that the pulsar is not able to produce the TeV 
gamma-ray emission si mply with an interstellar matter 
(ISJVI) density of 1 cm'^. iDubner et al. (2002) reported a 
denser ISIVI distribution of n ~ 10 cm-^, derived from the 
HI observations with ATCA, which is indeed valid only 
in the northwest radio limb. If such a high density was 
uniformly applicable, the total energy would be reduced 
to Wp = 3.2 X 10^° ergs. Dubner et al. (2002), however, 
also mentioned that the southeast radio limb of the SNR 
showed n ^ 0.4 cm"-^, which would yield to a higher Wp. 
In any case, it is revealed that the spin-down energy can 
not drive the acceleration of protons during its character- 
istic age, even if the uncertainty in the distance (see § 1) 
is taken into account. Therefore the hadronic gamma- 
ray production originated by the pulsar spin-down en- 
ergy was strongly unlikely. Note that bremsstrahlung 
emission is not dominant in the TeV gamma-ray band 
under these ambient densities. 

5.3. Leptonic scenario 

Secondly, we discuss the leptonic scenario, and here 
we used a simple one-zone IC model to reproduce the 
multi- wavelength spectra. 

Case I: The muffi-band SED is plotted in Fig. [H The 
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10"' 10"' 10"' 10^ 10^ 10 10^ 10° 10' 10' 10" 10" 10^° 

Energy [eV] 

Fig. 6. — SED and leptonic model curves. The references of data 
are given in Table [2] Assumed electron spectrum was a broken 
power-law and the derived model parameters are listed in Table |3] 
The IC spectra on each of CMB, IR, and star light are represented 
by dotted, dot-dashed and dashed curves, respectively. The pion 
decay model curve in Fig. \5\ is shown by the thick curve for com- 
parison. 

data points derived from this work are represented by 
filled triangles and references to others are summarized 
in Tab le [H Since BeppoSAX/PDS data in IJ VIineo et all 
(|2001t ) should be corrected by the intercalibration factor 
of about 20% between MECS and PDS, here we plotted 
the corrected PDS data (Mineo 2007, private communi- 
cation). The highest energy bin of the PDS represents 
a 1 cr upper limit. IBIS measurement reported extended 
emission with a possible spectral cutoff at about 160 keV, 
at the 2.9 a confidence level from the extrapolation of the 
power-law spectrum. Then we present an upper limit in 
Fig. [5] derived from multiplying the 1 a upper limit by 
a factor of 3. Note that, though the data analysis of 
the coded-mask instr ument of IBIS is ideally designed 
for point-like sources. [Porot et all (120061 ) have extracted 
the spectrum of the observed exte nded emission follow- 
ing a suited method developed bv^Rena ud et al.l (|2006f ). 
Arrows in the radio band show the whole emission from 
MSH 15-52 or RCW 89 listed in ldu Plessis et al.l (fTOQSh . 
which we treated as upper limits for the faint PWN. As 
well as in the case of the hadronic scenario, the COMP- 
TEL and EGRET fluxes were not used as they contain a 
component of pulsed emission (Kuiper ct al. 1999). The 
CMB field density on its own could not account for the 
TeV gamma-ray fiux, and hence we ad ded IR and opti- 
cal (s tarlight) photon fields, following lAharonian et al.l 
(|2005f ). Here we used an interstellar radiation field 
(IS RF) derived f rom the la test (v50p ) GALPROP pack- 
age (|Porter k StronEil2005t [Strong fc Moskalcnko 2006). 
The ISRF is given for three components (CMB, IR from 
dust, and optical starlight) as a function of a distance 
from the Galactic center R (in kpc) and distance from 
the Galactic plane z (in kpc). We extracted the spectra 
at (R,z) — (5.6,-0.11), which is shown in Fig. [7l and 
the nominal value is 1.4 eV/cc both for IR and optical 
light. The ISRF field densities did not change drasti- 
cally when we considered the ambiguity in the distance, 
within 1.4-0.89 eV/cc and 1.4-0.90 eV/cc for IR and op- 
tical, respectively. Note that these values extracted from 
GALPLOP do not represent the local densities. 
The best-fit curve is shown in Fig. [6] with the assump- 
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Fig. 7. — Interstellar radiation field from the GALPROP pack- 
age (v50p) at {R, z) = (5.6,-0.11) kpc. From lower frequencies, 
CMB (dotted), IR emission from interstellar dust (dot-dashed), 
and optical photons from stars (dashed) are shown. The solid line 
represents the sum of the three components. 

tion of a broken power-law electron spectrum; 



1 + {EjE^,,p2-'n) 



exp i~Ee/E 

n 



(5) 



where E^. is the electron energy, 71 is the spectral index of 
the injected electrons, 72 is that of the cooled electrons, 
and Ehi- and Smax are the break and maximum electron 
energies, respectively. The derived parameters are listed 
in Table H 

Fig. [8] shows a close-up of the SED and the de- 
rived model curves. According to the IBIS measure- 
ment, the unpulsed emission is dominated by the PWN 
(|Forot et al.l I2006D . which is also supported by the 
fact that the timing analysis showed unpulsed emission 
from the pulsar is indeed several factors le ss than the 
pul sed component (ICusumano et all 120011 ; iForot et "aLl 
12006:). The model curve well reproduced the hard X-ray 
spectra. The region closed with thick solid lines in Fig. [5] 
represents Chandra measurement excluding the pulsar, 
which corresponds to the sum of the " diffuse PWN", 
"jet" and "outer arc" fluxes listed in iDelanev et all 
(|2006r ). The FOV of the Chandra measurement was 
smaller than the MECS signal region, which may cause 
the apparent discrepancy in the soft X-ray band. There- 
fore we didn't take care of the data for the fit. 

The obtained magnetic field of 17 fiG was consis- 
tent with Ahar onian et al. (2 005), and higher than pre- 
vious indications ( Gaensler et all I200l . and references 
therein). The IR densit y was also compatible with that of 
lAharonian et all (|2005f ). Although the electron index of 
7i = 1.2 is much harder than the v alue of 2 pr edicted by 
gener al acceleration theories (e.g., iBlandford & EichleJ 
|1987( ). the corresponding radio spectral index of F^, oc 
i/"'' , ar = —0.1, was within typical value s for PWNe, 
-0.3 < < (jWeiler fc Panagial Il978^ . The ratio 
of observed X-ray to gamma-ray flux showed the syn- 
chrotron cooling was dominant compared to IC cooling. 
Thus, the spectral break was due to the synchrotron cool- 
ing, -Bbr.cooi — GTrrn^c^ / {B^ aTT) , where i? is a magnetic 
field of the PWN and (Tt is the Thomson cross section. If 
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Fig. 8. — Chandra flux is represented by a region closed with 
thick lines. Other data and the synchrotron and IC model curve 
in Fig. [6] are also shown. 
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Fig. 9. — SED and the model calculation of the leptonic scenario 
with a single power-law electron distribution. The representation 
of each component is the same as in Fig. [6] 

we accept a magnetic field of 17 /j,G, the break in the pho- 
ton spectrum was expected to occur at '^0.2 keV and ^ 
1 eV for the age of r = 1700 yr and 20kyr, respectively. 
However the electron break energy i?br of 77 GeV, de- 
rived from the model fit, predicts a break in the photon 
spectrum at lower energy of 2 x 10~^ eV, which sug- 
gests larger B or r. By integrating eq. ^ above 1 GeV, 
the total energy We amounts to 3.0 x lO''^ ergs. In the 
case of T = 1700 yr and tq = 30 yr, ~4% of i?tot is re- 
quired to be converted to — an age of 20kyr may 
more easily supply sufficient energy. 
Case II: The expected break in the photon spectrum at 
^ 0.2 keV, as mentioned above, indicates that the cooling 
process might not be effective. Then we could alterna- 
tively apply a single power-law electron spectrum using 
eq. ([5]) with 71 = 72. The reproduced SED is shown 
in Fig. [9] with parameters in Table [3l A radio spec- 
tral index — —0.6, derived from the electron index 
of 71 = 72 = 2.2 is close to a mean value for SNRs of 

ar 0.5 fe.g.. lGreenlll991[ ) rather than for PWNe of 

—0.3 < < 0. Here, the IR energy density is more 
than 3 times higher than that of GALPROP. There are 
several possible factors for variations in the local ISRF 



density of this magnitude, including uncertainties in the 
model itself but also on the conditions in the environment 
of MSH 15-52. IRAS 15099-5856 is a candidate local 
source for an increased IR field, if it is associated with 
PSR B1509— 58. IR from the s urrounding dust grains 
were also suggested by du Ple s sis et al.l ()l995ll . In Fig. [9] 
the sub-TeV gamma-ray flux was not reproduced. If we 
attempt to reproduce it, the IC peak energy, which was 
determined by the maximum electron energy, must lie 
below ~0.1 TeV. The effect of a stronger magnetic field 
would be to reduce the IC peaks to lower energies and 
also suppress the peak height. Since these IC spectra 
were dominated by the IR component, if the IR field 
density is increased, the IC peak is raised vertically in 
the SED. Therefore the combination of a strong mag- 
netic field and IR radiation field, such as B ~lmG and 
UiYi '-~^7000eV/cc, could reproduce the TeV gamma-ray 
spectra. However such a strong magnetic field would 
result in sufficient synchrotron cooling and cause a pho- 
ton spectral break at '^0.06eV, in conflict with the as- 
sumption of this single power-law model. Although the 
requirements from energetics are looser than that of the 
broken power-law model, the SED is not well reproduced. 

Finally, in both cases we can calculate the synchrotron 
and IC luminosities from the obtained model curve: 
16(case I)-12(case II)% and 0.6-0.4% of the current spin- 
down energy E(t) is radiated via synchrotron emission 
and IC, respectively. It is also possible to estimate the 
equipartition magnetic field strength, B^q^ from the fit 

results as i?oq — \/^nWe/V, where V is the volume 
of the emission region. For Case I, if we adopt V = 
1.3 X 10^8 cm3 in lSeward et all (fl98l . Beq = 76 jliG is ob- 
tained, which means the PW N is particle-dominated as 
suggested in Chevalier I ()2004f l . Assuming V as a sphere 
at 5.2 kpc with 17' radius (corresponding to the MECS 
signal region), B^q = 6/xG is obtained, which is compa- 
rable with the previous indications (e.g. iGaensler et al.l 
[2OOI ). In this case the PWN is poynting-flux domi- 
nate d. The Crab nebula is nearly in the equipartition 
(e.g.. I Atovanll 1 999t l . while several PWN e have a deviation 
with a few factor from equipartition (jChevalier 1 12005| ) 
and for other PWNe the situation is not clear at present. 
As for the energetics, the total energy we have derived, 
We of 3.0 - 0.5 X 10^8 ergs and Wp of 3.2 x lO^o-si ergs, 
depend on the distance and the age of the pulsar. The 
uncertainty in the distance would modify the total en- 
ergy by only a factor of a few, not an order of magnitude. 
Our estimations have been optimistic and the energetics 
would require a higher efficiency VFeZ-E'tot when expan- 



sion loss of PWNe is considered (IPacini &: Salvati 



Chevalier fc FranssonI fl993 Ivan der Swaluw et al 



1973 



2001 



Bejger fc H acnscl 2001 IChevalier 1120051 ). Although sim- 
ilar studies concerning the source of radiating parti- 
cles have recent ly been started for other PWNe (eg., 
iFunk et "al] l2007l la) . conclusive solutions have yet to be 
reached. 

6. CONCLUSION 

CANGAROO-III observed the SNR MSH 15-52 con- 
taining PSR B1509-58 for 48.4 hours in 2006 and de- 
tected VHE gamma-ray emission at the 7a level. The 
obtained differential flux and the intrinsic extent of TeV 
gamma-ray emission are consistent with the previous 
H.E.S.S. result. Studies of the multi-wavelength spec- 
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tra was performed, based on both hadronic and leptonic 
models. In the leptonic scenario, an IR photon field and 
a cooled broken power-law spectrum of electrons are nec- 
essary to reproduce the TeV gamma-ray emission. From 
the point of view of the energetics, if we do not take into 
account the expansion loss, a typical supernova could 
provide sufficient energy for electrons to reproduce the 
SED, while hardly to protons. The morphology of the 
TeV gamma-ray emission, however, does not support the 
supernova explosion as the global energy source. Elec- 
trons can also be accelerated enough to reproduce the 
SED when > 4% (tq = 30 yr is assumed) of the rota- 
tional energy is supplied to the kinetic energy. If the 
pulsar was older than its characteristic age of 1700 yr, 
e.g., 20kyr, the required acceleration efficiency would be 
reduced. 

Filling in the gaps in the SED is very important for 
the discussion above. IR/optical measurements of the 
PWN emission are crucial to determine the synchrotron 
spectrum and the electron scenario. Additionally, the 
determination of the spectral break in the synchrotron 
or IC component would help to resolve the long-standing 
question of the age of this complex system. We await 
the results from al l sky survey of the rec ently launched 
IR satellite Akari (iMurakami et al.ll2007l) . GLAST and 
the next generation of the ground-based lACTs such 



as CTA (iHermann et al.l[2007l) and AGIS ()AGIS|[2001 

are expected to determine the IC spectra. In addition, 
X-ray observations of the whole PWN, excluding PSR 
B1509— 58, are required in order to accurately estimate 
the flux of synchrotron nebula emission. 
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TABLE 2 

Summary of data used in the SED analysis. 



Detector 


Marker 


Pulse component 


Spatial component 


Reference 


ATCA 


open circle 


unpulsed 


N 


{1)(2) 


Chandra/AClS 


closed region(thick lines) 


pulsed+unpulsed 


NJ 


(2)(3) 


SAX/MECS 


closed region(thin lines) 


unpulsed 


PNJC 


(4) 


SAX/VT)S 


open square 


unpulsed 


PNJRCO 




INTEGRAL/lBlS 


filled circle 


unpulsed 


PNJRC 


(6) 


RXTE/PCA+HEXTE 


closed region(dotted lines) 


unpulsed 


PNJRCO 


(7) 


COMPTEL 


open diamond 


pulsed+unpulsed 


PNJRCO 


(8) 


EGRET 


open cross 


pulsed+unpulsed 


PNJRCO 


(8) 


H.E.S.S. 


open triangle 


pulsed+unpulsed 


PNJRC 


(9) 


CANGAROO-III 


filled triangle 


pulsed+unpulsed 


PNJRC 


this work 



References. — (1) Gaensler et al. 1999; (2) Gaensler et al. 2002; (3) Delaney et al. 2006; (4) Mineo et 
al. 2001; (5) Mineo 2007, private communication (6) Forot et al. 2006; (7) Marsden et al. 1997; (8) Kuiper 
et al. 1999; (9) Aharonian et al. 2005; 

" N: Diffuse PWN define d in ref (2), J: Jet an d outer arc defined in ref (2) and (3), P: Pulsar, C: Central 
Diffuse Nebula defined in lTrussoni et al.l 1 1199^) , R: RCW89, O: Outside MSH 15-52 



TABLE 3 

Summary' of parameters used in leptonic model. 



Electron B 71 72 Sbr -Emax Uin We Lsync'' Lic° 

spectrum [/^G] [GeV] [TeV] [eV/cc] [ergs] [ergs/s] [ergs/s] 

broken P.L. 17 1.2 2.7 77 2.5 X 10^ 2.3 3.0 X 10"'® 2.9 X 10^^ 1.2 X 10^^ 

single P.L. 20 2.2 2.2 - 1.3 X 10^ 4.5 5.4 X 10"^'^ 2.1 X 10^^ 7.7 X 10^"^ 
Luminosity of synchrotron emission calculated from the model fit curve. 
^ Luminosity of inverse Compton emission calculated from the model fit curve. 
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